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Abstract — THIS PAPER IS ELIGIBLE FOR THE BEST STU- 
DENT PAPER AWARD. While the efficiency of MIMO transmis- 
sions in a rich scattering environment has been demonstrated, less 
is known about the situation where the fading matrix coefficients 
come from a line-of-sight model. In this paper, we study in 
detail how this line-of-sight assumption affects the performance 
of distributed MIMO transmissions between far away clusters of 
nodes in a wireless network. Our analysis pertains to the study 
of a new class of random matrices. 

I. Introduction 

The aim of the present paper is to study the number of 
spatial degrees of freedom of multiple-input multiple-output 
(MIMO) transmissions in a wireless network with homoge- 
neously distributed nodes, under the following classical line- 
of-sight propagation model between node k and node j in the 
network: 

h 3 k = ■ (1) 

In the above equation, A is the carrier wavelength and rjk is 
the internode distance. From a mathematical point of view, 
these matrices are interesting objects, as they are halfway 
between purely random matrices with i.i.d. entries and fully 
deterministic matrices. Indeed, the internode distances rjk are 
random due to the random node positions, but there is a clear 
correlation between the matrix entries. 

Let us recall that the degrees of freedom of a MIMO trans- 
mission are defined as the number of independent streams of 
information that can be conveyed simultaneously and reliably 
over the channel at high SNR. Under the assumption of a 
channel fading matrix H with i.i.d. entries, this number of 
degrees of freedom is directly proportional to the number of 
antennas used for transmission and reception lfT3l . 

The performance of MIMO systems in line-of-sight environ- 
ment has been analyzed by various authors (see e.g. 0, JTT]) 
in the literature. Our intention here is to study this performance 
in the context of wireless networks, where large clusters of 
nodes are used as virtual multiple antenna arrays. In this 
case, MIMO transmissions may not benefit from all possible 
degrees of freedom; it was indeed observed in [4 1 that under 
the above propagation model (Q}, MIMO transmissions suffer 
from a spatial limitation; if A denotes the network area, n 
the number of nodes in the network (assumed to be uniformly 
distributed) and A the carrier wavelength, then the number of 



spatial degrees of freedom of any MIMO transmission in the 
network cannot exceecf] 

min (n, VH/Xj . (2) 

In case the network area A remains reasonably large, this does 
not prevent the possibility of transmissions with full degrees 
of freedom in the network. Yet, transmissions between clusters 
of nodes confined to smaller areas and moreover separated by 
long distances may suffer from even more spatial limitations. 
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Fig. 1. Two square clusters of area A separated by distance d. 

In 0, (HI, it was shown independently that for two clusters 
of area A separated by distance d, as illustrated on Fig. 1, 
at least the following spatial degrees of freedom could be 
achievec0: 

min (n, VA/Xj , when 1 < d < VA, 
min (n, A/Xdj , when VA <d< A/X. 

The situation is summarized on the graph below: 
dof 
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Fig. 2. Spatial degrees of freedom between two clusters 
of area A separated by distance d. 

'See (4) for a precise statement. 

2 See Section |n] Theorem ff] for a precise statement. 



We see that the lower bound ® matches the upper bound 
(O in the case where the inter-cluster distance is smaller than 
or equal to the cluster radius (d < yA), but nothing similar 
holds for d > \fA. Our aim in the present paper is to close 
this gap and to show that in the regime where d > \f~A, the 
actual spatial degrees of freedom of the MIMO transmission 
do not exceed those found in Q (up to logarithmic factors). 
As a corollary, this would imply that when d > A, the number 
of degrees of freedom is bounded by 1 . 

In order to show this, we rely on an approximation whose 
validity is not fully proven here; it is however discussed in 
detail at the end of the paper. Our approach leads to an 
interesting result on the asymptotic behavior of the spectrum 
of random matrices that appear not to have been previously 
studied in the mathematical literature. 

II. Spatial degrees of freedom 

Let us consider two square clusters of area A separated 
by a distance d, one containing n transmitters and the other 
containing n receivers uniformly distributed in their respective 
clusters, as illustrated on Fig. 1. We are interested in estimat- 
ing the number of spatial degrees of freedom of a MIMO 
transmission between the two clusters: 



Yj =Y,^Fh 3k X k + Zj, 3 = l,. 

k 



, n, 



where F is Friis' constant, the coefficients hj k are given by 
the line-of-sight fading model (Q3 and Zj represents additive 
white Gaussian noise at receiver j. The distance rj k between 
node j at the receiver side and node k at the transmitter side 
is given by 



r jk = ^]{d + \fA (xj + w k )) 2 + A (y 3 - z k f 



(4) 



where Xj ,w k ,yj, z k <E [0,1] are normalized horizontal and 
vertical coordinates, as illustrated on Fig. 3. 
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Fig. 3. Coordinate system. 

Assuming full channel state information and perfect coop- 
eration of the nodes on both sides, the maximum number of 
bits per second and per Hertz that can be transferred reliably 
from the transmit cluster to the receive cluster over this MIMO 
channel is given by the following expression: 



Cn 



max log det( / + HQ H*). 

Q>0:Q kk <P,\fk 



where Q is the covariance matrix of the input signal vector 
X = (Xi , . . . , X n ) and P is the power constraint at each node 
(in order to simplify notation, we choose units so that the other 
parameters, such as Friis' constant F, the noise power spectral 



density An. and the bandwidth W do not appear explicitly in 
the above capacity expression). 

In the sequel, we make the following two assumptions: 

1) d, A both increas^l with n and satisfy the relation 
VA < d < A/X, which is the regime of interest to us (see 
Fig. 2). 

2) P = {d+\/~A) . b ecause me average distance between 



two nodes in opposite clusters is d + v A and because the 
MIMO power gain is of order n, this power constraint ensures 
that the SNR of the incoming signal at each receiving node 
is of order 1 on average, so that the MIMO transmission 
operates at full power. Imposing this power constraint allows 
us to focus our attention on the spatial degrees of freedom of 
the system. 

By choosing to transmit i.i.d. signals (i.e. taking Q = PI), 
we obtain 

C n > logdet(J + PHH*) 

and using Paley-Zygmund's inequality, the following result 
was further shown in 

Theorem 1. Under assumptions 1) and 2), there exists a 
constant K\ > such that 



C n > logdet(7 + PHH*) > K x min 



with high probability as n gets large. 



A/Xd \ 
\og(A/Xd)J 



This result shows that the number of spatial degrees of 
freedom of the MIMO transmission can reach A/Xd (up to 
a logarithmic factor), when the number of nodes participating 
to the MIMO transmission is large enough. 

As mentioned in the introduction, a natural question is 
whether it is possible to find a corresponding matching upper 
bound on the capacity. In order to answer this question, let us 
first observe that any matrix Q satisfying the above constraints 
also satisfies Q < nPL Thus, 

n 

C n < log det(/ + nPHH*) = ^log(l + A fc ) (5) 

k=\ 

where Ai > A2 > . . . > A„ are the eigenvalues of nPHH*. 
The number of significant eigenvalues of nPHH* therefore 
determines the number of spatial degrees of freedom. The 
direct analysis of these eigenvalues appears to be difficult, so 
we proceed by approximating the matrix nPHH* by another 
matrix GG* , easier to analyze. 

Claim 2. Let m = A/Xd and G be the matrix whose entries 
are given by 

(6) 



— p-ZirimyjZk 

yjk — c , 



3 By "increasing with n", we mean that A 
powers 0, 7 > 0. 



n,P and d = n 1 for some 



where yj, Zk, 1 < j, k < n are the same random variables 
as in expression (0). Then under assumptions 1) and 2), the 
following approximation holds: 

logdct(/ + nPHH*) = logdet(7 + GG*) (1 + o(l)) 

with high probability as n gets large. 

We discuss this approximation in detail in Section [HI] For 
the time being, observe first that by expression (0, the above 
approximation is equivalent to saying that the number of 
significant eigenvalues of nPHH* and GG* do not differ 
in order as n gets large. Some numerical evidence of this fact 
is provided on Fig. 4 for a given set of parameters (a similar 
behavior is observed for a wide range of parameters). 
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Fig. 4. Eigenvalues of nPHH* (blue) and GG* (red) for the 
parameters n = 500, A = 10'OOOto 2 , d = 300m, A = 0.1m 
(so to = A/Xd ~ 333). 

It can be observed on the figure that the eigenvalues drop 
to zero after a threshold of order m — A/Xd for both matrices 
nPHH* and GG* . Notice that because of assumption 1), we 
have m — n s , where the power 8 is a number greater than 0. 

The rest of the present section is devoted to the proof of 
the following statement. 

Theorem 3. Let to = A/Xd be such tha^ to ^> *Jn. Then 
under assumptions 1) and 2), there exists a constant K% > 
such that 

logdet(/ + GG*) < K2 min(n, to) logn 
with high probability as n gets large. 

This result shows that the lower bound found in Theorem 
Q] is tight (provided Claim [2] holds true and to 3> \/n), which 
is saying that the number of spatial degrees of freedom of a 
MIMO transmission between two clusters of area A separated 
by distance d is of order to = A/Xd, up to logarithmic factors. 

Let us also mention that applying the same technique as 
in |H0, the same matching lower bound on logdet(7 + GG*) 
can be found. This result on matrices G of the form © is 
interesting in itself, as these do not appear to have been studied 
before in the random matrix literature. 

Proof: First, observe that applying the same method as 
in JS] Lemma 2.2], the following concentration result can be 

4 i.e. m = n s , where 8 > 1/2. 



shown: for all e > 0, there exists some constant K > such 
that 

I logdet(7 + GG*) - E(logdet(7 + GG*))\ < Kn^ 2+e 

with high probability as n gets large. As to ^ %fn by 
assumption, what remains to be shown is that there exists a 
constant K2 > such that 

E(logdet(J + GG*)) < K 2 min(n, to) logn 

as n gets large. Observe that this was the only part of the proof 
that requires m \pn. It follows that a sharper concentration 
bound would immediately yield a stronger result in Theorem 

El 

In order to upperbound E(logdct(7 + GG*)), let us now 
expand the determinant: 

E(log det(7 + GG*)) 

log U + £ det(Gj xn G* Jxn ) 

^ fc=l JC(1 n} 

J7| = fe 

/ " 

< log 1 + £ ]T E(det(G Jxn G* Jxn )) 

^ fe=l JC(1 71} 

1.7!=* 

where we used Jensen's inequality. Using the fact that the yj 
are i.i.d., we further obtain that E(det(G i 7 X nG^ xrt )) only 
depends on the size k of the subset J, so 

E(logdet(7 + GG*)) 



< log (l + Q E ( det GkxnG* kxn 



log ( 



k=l 



E 



^ det(G fex xG£ XI ) 



IC{1 0} 

\I\=k 



= log 1 + E(det(G fcxfc G^ xfc )) 

^ k=l ^ ' 

where we have used this time the Cauchy-Binet formula 
together with the fact that the are i.i.d. We thus see that 
in order to upperbound E(logdct(7 + GG*)), it is enough 
to control E(det(Gfc X fcG£ xfe )), where Gkxk is the upper left 
k x k submatrix of G. 

We will show that, similarly to what has been observed 
numerically for the eigenvalues Xk, E(det(G/ cx / c G^. xfc )) drops 
rapidly for k greater than a given threshold of order m, which 
will imply the result. 

Using the definition of the determinant, we obtain 

E(det(G fexfe G£ xfe )) 



aes k 



which in turn leads to 

E(det(G kxk G* kxk )) 



a£S k 



3 = 1 

k 



k 'S} ^ ^JJ 27rim ( z j - Mi)) 

f f ( \ — p—2mm(zj — Zi) ' 

k\E z \ dot 



E z det 



2wim(z j -z l ) J!<^< fc 

Multiplying row j by e 7rimZj and column / by e^ 7 "™ 2 ^ we 
reduce the problem to computing the following determinant: 

sin(7rm(zj — Z;)) 

nm{ Zj - zi) ) 1 < j>l < ky r 

Operators and Fredholm Theory. The key observation is that 
the above expected value of the determinant can be seen 
as a classically studied quantity in the Fredholm theory of 
integral operators. This allows us to deduce precise estimates. 
A reference for the material discussed below is ||9"1 . 

Consider the continuous function K m (x,y) = ""j^ff 
on [0, l] 2 and the associated operator K m : C([0, 1]) — > 
C([0,1]) defined as 

,/ x f 1 sm(m(x — y)) , , . , 



o tt(z - y) 

The p th iterated kernel K p of an operator K is defined as 

K 1 = K and 

K p (x,y)= f K p - 1 (x,z)K(z,y)dz 
Jo 

Associated to this is the p th trace of K: 



.4 



p 



K p (x,x)dx 



Define as well the compound kernel Ky\ e C([0, l] 2p ) as 

(K{xi,yi) K(x 1 ,y 2 ) ■■■ K(x 1 ,y p y 
K[ p ] (x, y) = det : : ■ ■ . : 

\K(x p ,yi) K(x p ,y 2 ) ... K(x p ,y p ) J 

for x = (x%, . . . , x p ) and y = (yi, . . . , y p ). In this notation, 
the quantity we are interested in is 



E z ^det 
where 



nmizj - zi) J x <^< k ) J m k 



d k = 



[QA] k 



K m< [ k ](x,x)dxi ■■■dx k . 



Since K m is a compact operator, it has a discrete spectrum 
A*i > 1^2 > ... It is immediate to see that A p — ^ /x p . 
Furthermore, the quantities d k and A p are related by the 
following recurrence relation, which follows from expanding 



the determinant in the definition of K[ p ] and regrouping equal 
terms (see [9]). 

k 

kd k = ]T(-1)p- 1 A p 4- p . 
p=l 

Using the fact that A p = Vi, it can be seen that the only 
solution to the above recurrence (with do = 1) is 



d k 



E 



7 , Mil /''.- ' ' ' /''. • 



(7) 



We conclude that it is sufficient to estimate the eigenvalues of 
the operator K m in order to upperbound d k . 

Since the kernel K m is translation invariant, it can be 
defined equivalently on [—1/2, 1/2], and this new operator has 
the same eigenvalues. This operator is called the sine kernel 
and is well known in signal processing, since it is the Fourier 
transform of the indicator function. It was originally studied by 
D. Slepian in [12| and precise estimates exist on the behavior 
of its eigenvalues. In particular, we have the following recent 
result from [1. Theorem 3]: 

Theorem 4. Let 6 > 0. There exists M > 1 and c > such 
that, for all m > and k > max(M, cm), 

fl k < e -<5(fc-cm). 

This theorem essentially says that the eigenvalues [i k decay 
exponentially for k > cm. The direct consequence of this is 
that d k decays like exp(— 6 (k — cm) 2 /2) for k > cm, as we 
show below. 

Indeed, it follows that if we take k sufficiently large 
(i.e. larger than cm), the sum in © always contains at least 
one term with exponential decay. Recalling that all eigenvalues 
\x k are bounded by 1, we obtain 

d k = " ' 

ii<...<i k 

< min(e- 5(ll - cm \l)---min(e-^- 1 - cm \l) 

it <...<i k -i 

oo 

X e -S(i k -crn) ^ 

as i k is necessarily greater than or equal to k in the above 
sum. Now, define K = k — cm and observe that 



x p — S(k—cm) 

E e — 5{i k — cm) _ v 
1 - e~ s 

i k =k 



< Ce 



-8K 



for some constant C > 0. So by induction, we further obtain 

d k <C K f[c- s ^<C K e- 6K ^ 2 . 

This gives us the estimate we were after for d k : 

1, if k < cm, 



d k < 



C k e -S(k-cm) 2 /2 if k> cm. 



Gathering all estimates together, we finally obtain 

E(log det(7 + GG*)) 

(cm n \ 

1 + J2n 2k + J2 n 2k C k e-^ k ' cm ^l 2 \ 
k=l k=cm+l ' 

< (cm + 1) logn + 0(1), 
which concludes the proof of Theorem [3] ■ 

III. Discussion 

Our aim in the following is to provide a justification for 
Claim |2] Let us first recall the definition of both H and G: 

27rir jfc /A 

h jk = and g jk = e- 2 * im v^, 

where m — A/Xd and 

r jk = \J (d + VA( Xj + w k )) 2 + A (y 3 - z k ) 2 . 

Given the chosen power constraint P and the fact that d > 
\f~A, it follows that the amplitude of the normalized fading 
coefficent \JnP hj k is of order 1, matching that of gj k . Let 
us now compare the phases of these two coefficients. Using 
a Taylor approximation to quadratic order around in the 
variable (j/j — z k ), we get 

r jk ~d + s/A( Xj + w k ) + (A/d) ( yj - z k ) 2 /2 

= d+VA( Xj + w k ) + (A/d) (y]/2 + zl/2) - (A/d) y 3 z k . 

Hence, 

2nir jk /\ ^ T ._ 2iri(uj +v k -(A/ \d) yjZ k ) 

e — IL J k • — e ; 

where 

Uj = (d/2 + y/Axj + (A/d) y]/2)/\ 
v k = (d/2 + VAw k + (A/d) zl/2) I \. 

Notice moreover that the eigenvalues of HH* do not depend 
on the particular values of the Uj's or w^'s; they are therefore 
the same as the eigenvalues of GG* , which shows finally that 

VnP h jk ~ g jk , Vj,k. 

This entry-by-entry approximation adds some plausibility to 
Claim E] 

IV. Conclusion and perspectives 

The goal of this work is to give precise estimates on the 
number of spatial degrees of freedom in large MIMO systems 
in a line-of-sight environment. An upper bound for a model 
closely related to the line-of-sight model has been given, 
and the similarity of the models is supported numerically. 
As such, it remains to be shown that the eigenvalues of 
the two models are indeed very close, in order to bound 
| logdet(7 + nP HH*) - logdet(7 + GG*)\; this is work in 
progress. 



As a by-product, the spectral properties of random matrices 
G of the form gj k = e~ 2vunyjZk have been studied in this 
paper. These matrices are not unrelated to both Vandermonde 
matrices and random DFT matrices that appear in other 
contexts in the literature on wireless communications [7 1, [ 10|, 
03), lfl5l and compressed sensing (2), Gl, respectively: 

- Vandermonde matrices are simply obtained by choosing 
i/j = j /n deterministic instead of uniform and i.i.d. on [0, 1]. 
In this case, the matrix entries are given by e - 27rim i 2 fc/™. 
Our analysis technique does not allow us to reach the same 
conclusions directly for these matrices, but we believe this 
method to be fruitful for further developments. 

- Random DFT matrices are obtained from the above Van- 
dermonde matrices by further replacing mz\ , . . . , mz n by a 
subset {li, . . . , l m } of m integers chosen uniformly at random 
in {1, ... , n}, so the matrix entries become e - 27 ™j'*:/™. Again, 
our technique can be applied in this setup. 
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